Quantum Scattering of Distinguishable Bosons using an Ultracold Atom Collider 
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We describe a new implementation of magnetic collider for investigating cold collisions between ultracold 
atomic clouds in different spin states, and we use this to investigate scattering involving both even and odd order 
partial waves. Our method relies on the axial assymetry of a double- well magnetic trap to selectively prepare the 
spin state in each cloud. We measure the energy dependence of s, p and d partial wave phase shifts in collisions 
up to 300 pK between 87 Rb atoms in the 55 , i/ 2 , F = 1, ttif = —1 and 55 f i/ 2 , F = 2, vtlf — 1 states. 

PACS numbers: 34.50.-s, 03.65.Nk, 34.10,+x, 32.80.Pj 



Collisions in ultracold and degenerate quantum gases play 
a key role in many of their interesting properties [1]. So 
far, investigations with ultracold atoms have been mostly con- 
cerned with s-wave scattering processes, but now nonzero par- 
tial waves play a critical role in many investigations, (see, e.g., 
101). A magnetic collider scheme for determining the contri- 
bution made by higher-order partial waves was recently im- 
plemented lyl 01 • In these experiments the atoms were in the 
same spin state, limiting the collisions to those involving only 
even-order partial waves — a consequence of the particles be- 
ing indistinguishable bosonic particles. 

In the present work, we extend our collider method to dis- 
tinguishable bosons for which the scattering is fundamentally 
different since both odd and even angular momentum compo- 
nents are allowed. As in our original work [3], spin-polarized 
87 Rb atoms are loaded into a magnetic double-well poten- 
tial which is then transformed to a single well to initiate a 
collision. Here, however, one of the clouds is converted to 
a different spin state prior to collision making the scattering 



patterns crucially different. We observe the interference of 
s, p and d partial waves for collisions between atoms in the 
F = l,?7iF = — 1 and F = 2,m^ = 1 hyperfine ground 
states. Despite the complexity of the three- wave interference, 
we successfully determine the three partial wave phase shifts 
for energies up to 300 pK as measured in units of the Boltz- 
mann constant ks- 



The angular dependence of the two-body scattering prob- 
lem is described by the complex scattering amplitude f{6) 
101. Using the partial wave expansion, this is expressed as 
f(0) = ^ YZ^ + l )^ m ~ l)iMcos0), where P £ is 
the I th order Legendre polynomial and rje are the partial wave 
phase shifts which depend on the scattering potential and rel- 
ative wave vector k of the colliding atom pair. For the range 
of energies we focus on here, only the first three partial waves 
t = 0, 1, 2 contribute Q]. In this case the differential cross- 
section dcr/dVt = \f(6)\ 2 is given by 



J 



= —{sin 2 t?o + 9 sin 2 771 cos 2 6 H — - sin 2 772 (3 cos 2 — l) 2 + 6 sin 770 sin 771 cos(t7o — 771) cos 6 
ail k l 4 

+ 5 sin 770 sin 772 cos(t7o — 772) (3 cos 2 6 — 1) + 15 sin 771 sin 772 cos(t7i — 772) (3 cos 2 6 — 1) cos 6}. (1) 
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Because of the orthogonality and completeness of the Legen- 
dre polynomials, a fit of an interference expression in the form 
Eq. ([T} to a measured angular distribution directly gives the 
partial wave phase shifts 770, 771 and 772 irrespective of knowl- 
edge about absolute quantities such as particle flux 0] . 

Our experimental procedure is as follows. 87 Rb atoms in 
the 5Si/ 2 F = l,rriF = —1 (= |1)) state are loaded into 
a magnetic quadrupole-Ioffe-configuration (QUIC) trap 0] 
with trap frequencies u z /2tt = 11 Hz axially and uo p /2it = 
90 Hz radially. The details of loading the double- well trap and 
initiating a collision are much the same as described in 01. In 
summary, after rf-induced evaporation of the atoms to a tem- 



perature of approximately 2 pK we adiabatically transform 
the potential to a double well by raising a potential barrier 
along the axial dimension of the trap to split the cloud in half 
@]. The clouds are then further evaporatively cooled to a tem- 
perature of typically a few hundred nano-Kelvin, just above 
the Bose-Einstein condensation transition temperature. A col- 
lision between the clouds is initiated by rapidly transforming 
the potential back to a single well. The collision energy is 
selected by adjusting the well spacing in the double- well trap. 

To enable a collision between atoms in different spin states, 
we apply a two-photon pulse consisting of a microwave (~ 
6.8 GHz) and an rf (~ 2 MHz) photon (depending on the Zee- 
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FIG. 1 : (color online) After a collision between two atomic clouds in 
different spin states, pairs of diametrically-opposite scattered parti- 
cles will be distributed over the expanding Newton sphere according 
to the differential cross-section. Using a light beam resonant with 
only one of the states (depicted as red), an absorption image of the 
contribution of this particular state to the scattering halo is obtained. 




man splitting) to transfer |1) state atoms to the 5Si/2, F = 
2,mf = 1 (= 1 2)) state [10]. Due to the intrinsic axial asym- 
metry of the QUIC trap the clouds are situated at slightly dif- 
ferent magnetic field values immediately after the double- to 
single- well trap transformation. This enables us to selectively 
address and convert up to 90 % of the atoms in one of the 
clouds, while only 10 % of the atoms in the other cloud are 
converted to the |2) state. To first order, the |1) and |2) states 
have the same magnetic moment and experience the same con- 
finement potential. 

To selectively probe the scattered |2) state atoms we ap- 
ply a 20 [is pulse of resonant light on the 5Si/2, F = 2 — > 
5P 3 / 2 , F' = 3 transition along a radial direction shortly af- 
ter the end of the collision, and acquire an absorption image. 
This leaves the |1) state atoms undetected. An illustration 
of this is shown in Fig. [I] Alternatively, we can simultane- 
ously probe both the 1 1) and |2) state atoms by applying some 
5$i/2? F = 1 — > 5P 3 / 2 , F' = 2 light to pump all of the atoms 
to the F = 2 level shortly before the probing pulse. 

Figure |2| shows absorption images after a collision at 
E/ks = 135 /iK between atomic clouds in the |1) and |2) 
states. In Fig. |2a) only atoms in the |2) state have been 
probed, whereas in (b) atoms in both the |1) and |2) states 
are imaged. The distinct left-right asymmetry of the scat- 
tered atoms in (a) is the result of partial-wave interference 
between the odd (£ = 1) /7-wave and even s- and d- waves. 
The scattering amplitude of the /?-wave component changes 
sign at 6 = ±tt/2 as can be seen in Fig. [3] For the collision 



FIG. 2: (color online) Scattering images for a collision at E/ks = 
135 /iK, (a) probing only atoms in the F = 2 state and (b) probing 
both the F = 1 and F — 2 states. The asymmetry in the scattering 
pattern of (a) is due to /?-wave scattering. The corresponding angular 
scattering probability is shown in (c) with a fit to Eq. Q (solid line). 



energy in this example, where the J- wave contribution is rela- 
tively small, the /?-wave interferes constructively with s-wave 
for angles \6\ < tt/2 and destructively for |6^| > tt/2 where 6 
is defined with respect to the collision axis in the initial direc- 
tion of travel (i.e., for the |2) state shown in Fig.|21 \0\ < tt/2 
is to the left of the image). Since 6 is defined with the oppo- 
site sense for the |1) and |2) states, f(6) for the |1) state is 
complementary to that of |2) and imaging both states together 
results in a symmetric scattering pattern [Fig.|2tb)]. 

We analyze the absorption images of the scattering patterns 
using the method described in full . Briefly, we reconstruct 
the 3D distribution of the scattered atoms using the inverse 
Abel transformation fl2ll . The Abel-inverted image is divided 
into 30 angular bins which reflect the trajectories of scattered 
atoms in the harmonic potential. The number of scattered 
particles in each of the bins yields a measure of the angular 
scattering probability, which is proportional to the differential 
cross-section in Eq. |l}. We fit Eq. {I} to this data to obtain 
the partial wave phase shifts 770, Vi an d V2 for the s, p, and 
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FIG. 3: (color online) A graphical representation of the contributions 
to the scattering amplitude for the first three partial waves; (a) s- 
wave, (b) p-wave and (c) d-wave. The sign and magnitude of each 
£ term in f(6) is plotted vertically for a spherical scattering shell 
with a Gaussian profile in the radial direction. The relative scale of 
each is determined by rj£. In contrast to s- and d- wave, the /?-wave 
contribution to f(0) is antisymmetric in 0. 



^partial waves respectively. As emphasized by Buggle et al. 
0], this is an interferometric method which does not rely on 
absolute particle numbers and identifies only the amplitudes 
and relative signs of the phase shifts. The s-wave scattering 
length is known to be positive (repulsive interaction) for the 
states considered here so we choose the corresponding solu- 
tion where t]q < for our energy range. The collision energy 



E = mv* el /A = hk 2 /m is measured within a typical uncer- 
tainty of 5 fiK by determining the relative velocity v re \ from 
a linear fit to the position of the clouds over approximately 
2 ms either side of collision. In Fig.0]each phase shift value 
is the average of up to 10 measurements at the particular col- 
lision energy. The error bars on the data combine statistical 
uncertainty and errors associated with the fit to Eq. 

A comparison of the measurements to theoretical predic- 
tions is shown in Fig.|4| These are standard coupled-channels 
numerical calculations fl5ll for the collision of two 

atoms in hyperfine states F, M and F', M' in a low mag- 
netic field B with relative (partial wave) angular momentum 
£ and projection m. All channels {FM, F'M 1 £m} coupled 
by terms in the molecular Hamiltonian are included. Only 
channels with M tot — M + M' + m can couple to one 
another and because the collisions are from a single direc- 
tion (defined by the vector connecting the two initial sepa- 
rated atomic clouds) we need only include the M tot = 
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FIG. 4: (color online) The partial wave phase shifts for collisions 
between the |1) and |2) states. The symbols represent the s (A), p 
(O) and d (□) phase shifts extracted from the data and the solid lines 
are a theoretical calculation from a coupled-channels model. 



channels. The Hamiltonian contains the radial Tr and rota- 
tional T ro t kinetic energy terms, the electron-electron spin- 
spin interaction a 2 H ss (where a is the fine structure constant), 
the electron-nuclear spin interaction terms a 2 HM that gives 
the atomic hyperfine energies, and the strong chemical in- 
teractions described by the two adiabatic Born-Oppenheimer 
potential curves that correlate with two separated 2 S atoms. 
These potential curves correspond to the electronic states of 
1 S+ and 3 E+ symmetry. There are 8 s-wave channels needed 
to describe M tot = s-wave collisions of |2) and |1) atoms 
(5 open and 3 closed). There are also 18 M tot = p-wave 
collision channels (11 open and 7 closed), and 30 M tot = d- 
wave channels (18 open and 12 closed). All of these channels 
are included in the basis set for each partial wave. If the chan- 
nels are designated by the index j, so that the wavefunction 
for atoms in the entrance channel i is = ^ . \j)fji(R)/R 9 
the coupled Schrodinger equations, in a basis set defined 
by the separated atom quantum numbers, takes on the form 



h 2 d 2 f kl 

2fi dR 2 



E — E k 



fi 2 4(4 + i) 



fki(R) 



2fiR 2 

J2v kj (R)f j i(R)=0. (2) 



Here E k and Ik are the respective Zeeman energy and relative 
angular momentum quantum number of the two colliding sep- 
arated atoms for the magnetic field B, and the potential matrix 
elements Vkj define the interchannel coupling. These equa- 
tions are solved numerically using standard algorithms fl6ll . 
For comparison with the data, the calculation uses a magnetic 
field of 0.23 mT, and the scattering potentials are character- 
ized by a dispersion coefficient Cq = 4703 au and triplet 
a t = +98.96 ao and singlet a s = +90.1 ao scattering lengths 
consistent with O (1 au = E h a% where E h = 4.36 x 10" 18 
J and a = 0.0529 nm). 
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As can be seen in Fig.|4l our experimental observations are 
described well by the theoretical model. The dramatic change 
of the d- wave phase shift is a signature for the d- wave shape 
resonance known to occur for collisions between two 87 Rb 
atoms 00,0]. We estimate the position of the resonance 
to be (235 ± 50) /jK with a width of approximately 120 fiK 
(FWHM) from a Lorentzian fit to the data around the reso- 
nance. Calculated inelastic collision rate constants remain be- 
low 10 -13 cm 3 /s over the collision range of interest (com- 
pared to a maximum total elastic scattering cross-section of 
~ 1.6xl0 -11 cm 2 ), even when enhanced by the d-wave shape 
resonance. This is due to the exceptional case that both poten- 
tials have similar scattering phase shifts at low collision ener- 
gies for threshold 87 Rb spin-exchange relaxation Emm. 
Correspondingly, we do not observe any atom loss from the 
trap resulting from the collision. 

Two effects are not included in our analysis: state impu- 
rities in the clouds and the possibility of multiple scattering. 
The first of these is a difficult technical issue relating to our 
set-up and the second is of a more fundamental nature. With 
state impurities in both clouds, the collision processes which 
can occur are |1) + |2), |2) + |2), |1) + |1), and |2) + |1), with 
relative amounts depending on the density of impurities. If 
these effects were significant one would expect the presence of 
collisions between the |1) and |2) states in the "wrong" direc- 
tion to diminish the measured /7-wave contribution, whereas 
scattering due to the |2) + |2) and |1) + |1) collision pro- 



cesses would increase the perceived s- and d- wave phase shifts 
measured which is clearly not the case in Fig. |4] As for the 
second issue, we observe only approximately one-third of the 
total number of atoms scattered after a collision near the reso- 
nance, indicating that the probability of a secondary collision 
is relatively small. A detailed theoretical analysis of multiple 
scattering is difficult outside the s-wave regime, and partic- 
ularly near a d-wave shape resonance, since the energy and 
centre-of-mass of a subsequent collision depend crucially on 
the scattering angle after the first collision. 

In conclusion, we have investigated the energy dependence 
of collisions between two 87 Rb clouds in different spin states. 
Our experimental observations agree well with predictions 
from a theoretical coupled-channels model. We note that the 
collision between two such particles of different spins pro- 
vides a mechanism for producing spin entanglement. The re- 
sulting pair correlation could potentially be observed as in re- 
cent experiments on dissociating molecules ll22ll and colliding 
Bose-Einstein condensates Furthermore, the occurrence 
of a J- wave resonance and the resulting directionality of scat- 
tered particles may serve as a vehicle for the production of 
pair correlated beams. 
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